KPZ scaling in topological mixing 
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In the spirit of recent works on topological chaos generated by sequential rotation of infinitely thin 
stirrers placed in a viscous liquid, we consider the statistical properties of braiding exponent which 
quantitatively characterizes the chaotic behavior of advected particles in two-dimensional flows. We 
pay a special attention to the random stirring protocol and study the time-dependent behavior of 
the variance of the braiding exponent. We show that this behavior belongs to the Kardar-Parisi- 
Zhang universality class typical for models of nonstationary growth. Using the matrix (Magnus) 
representation of the braid group generators, we relate the random stirring protocol with the growth 
of random heap generated by a ballistic deposition. 

PACS numbers: 



Finding the ways of efficient mixing in liquids is one 
of the very popular subjects in physics of fluids because 
of numerous important technological applications rang- 
ing from the industrial production of blends and alloys 
to the percolation of turbulent flows of viscous liquids 
through porous media. In connection with that the ex- 
perimental and theoretical study of the chaotic motion 
of passively advected particles in two-dimensional flows 
in presence of moving rods ("stirrers") become very pop- 
ular during the last decade after the key work [l[ gener- 
alized and developed later in The most attention 
in these works is paid to the following question: how 
should the stirrers (playing the role of topological obsta- 
cles in the liquid) move to produce the strong chaos in 
passively advected particles in viscous liquid. This ques- 
tion is known as finding the best "stirring protocol" . The 
presence of the strong chaotic behavior due to motion of 
stirrers is guaranteed by topolo gica l arguments based on 
the Thurston-Nielsen theory fol, 1 1 dl] which postulates the 
existence of regions in the flow with pseudo-Anosov dy- 
namics when every point of the flow demonstrates the 
exponential stretching with positive Lyapunov exponent. 

On the basis of topological arguments it has been 
pointed out in [l| that the chaotic behavior in the flow 
of advected particles can be quantitatively characterized 
by measuring the Lyapunov exponents of the product 
of matrices representing the sequential permutation of 
neighboring stirrers in space. To be more precise, the fol- 
lowing "discretized" model has been considered (see, for 
example, [ill ? ]). Suppose that the viscous liquid passes 
through the comb of stirrers and in the same short inter- 
vals of time two neighboring in the comb stirrers permute 
carrying along some part of a liquid including a trace of 
advected particle. In the most transparent way such a 
model has been formulated in Q. In this work the liq- 
uid passed through the comb of n aligned infinitely thin 
rods. The permutations of neighboring rods leads to the 



entanglement of the trace of advected particle and can 
be characterized by a generator of the braid group B n+ i 
(the definition is given below). Thus, the time-ordered 

product of T sequential permutations, Wt ='Y\t=i °«t : ^ s 
a word written in terms of elementary operations - per- 
mutations of neighboring stirrers. Each generator of the 
braid group, o~\ (1 < i < n) representing the permuta- 
tion of stirrers i and i + 1 in the comb can be written 
in the matrix (Burau) representation (see, for example, 
[HI]). Thus, the word Wt is an n x n square matrix Wt 
containing the full information about the "stirring pro- 
tocol". Computing the Lyapunov exponents A{Wr} for 
different sequences of permutations and looking at the 
dependence A(T) one can conclude about the properties 
of the mixing regime. 

In this letter we formulate the mixing problem in a 
slightly different manner paying the most attention to 
its statistical aspects. Briefly, the outline of our consid- 
eration is as follows. We suppose that the matrix Wt 
is the time-ordered product of T generators of the braid 
group taken in a random order with the uniform distri- 
bution over the set generators {cti,CT2, ...,cr n }. In terms 
of the works P, 0, Q this procedure of mixing can be 
regarded as a random stirring protocol. For each gen- 
erated sequence Wt (for fixed T) we compute the asso- 
ciated Lyapunov exponent, A{Wt}- Averaging A{Wr} 
over the ensemble of all possible sequences Wt we find 
the mean value, A(T); also we compute the variance, 

u 2 (T) — (\{Wt} — A(T)) 2 . Our main result consists 
in the following. We show that in the broad interval 
of values T and n, the variance u(T,n) demonstrates 
the Kardar-Parisi-Zhang (KPZ) scaling (see, for review 
0), i.e. u(T,n) = n^ 2 g(T/n 3 / 2 ) with g(u) ~ it 1 / 3 for 
u <C 1 and g{u) ~ const for b»1. Such a behavior is 
typical for many non-stationary growth problems. It is 
believed now that KPZ universality class for correlated 
stochastic processes is as typical as the Gaussian statis- 
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tics for uncorrelated ones. 

To make the content of this letter as selfconsistent as 
possible, it is instructive to introduce the basic defini- 
tions and precisely formulate the model under consider- 
ation. We begin with the definition of the braid group, 
B n . The group B n+ \ has n generators <j\,...,<j n (and 
their inverses) with the commutation relations: 



(j i o i +\o 1 = a l+ i(T l (T i+ i 1 < i < n 



k\>2 



(1) 



Any arbitrary word written in terms of letters-generators 
from the set {cti, ...cr n , erf 1 , cr" 1 } — gives a particular 
braid. The graphic representation of the braid generators 
is shown in Figflji. The length of the braid is the total 
number of used letters. Diagrammatically the braid can 
be represented by a set of crossing strings going from 
the top to the bottom after subsequent gluing the braid 
generators. 



1 i-l i i+l i+2 n 

X 

1 i-l i i+l i+2 n 

Figure 1: Braid group generator en and cr 1 . 



The generators of the groups B n admit the matrix rep- 
resentations. For our purposes it is more convenient to 
use for the braid group generator Oi the Marnus repre- 
sentation instead of the Bureau one [l2|]. Thus, we have 
for a*: 
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<— row i: A = 




(2) 

For different i's (2 < i < n - 1 ) the 3 x 3 blocks A in © 
" slides" along the diagonal, while for i = 1 and i = n the 
corner blocks arc correspondingly A' and A" , where 



A' = 



-u 1 
1 



A" = 







(3) 



Note, that the matrices A' and A" are the generators of 
the simplest nontrivial group B% . It is known (l2l | that B3 
for u = — 1 coincides with the modular group PSL(2, Z), 
i.e. the group B3 is the central extension of PSL(2, Z). 
To speak rigorously, from this point of view, the matrix 
representation of -B3 used in |l| and in the consecutive 
works deals with the group PSL(2, Z), but not with B3. 



Now we are in position to precisely formulate the model 
under consideration and the main question of our in- 
terest. In what follows we consider for simplicity the 
semigroup B^ +1 , i.e. we construct the words of letters 
(generators) from the set {o~±, , , ,<J n } only (without the 
inverses). The extension of our consideration on the full 
group and related difficulties are discussed briefly at the 
end of this letter. 

Define the random walk on a semigroup B^ +1 . Let 
fx = — be the uniform measure on the sets of generators 
of B^ +1 . Consider the (right-hand side) random walk on 
Bn+i w ith fixed framings, i.e. regard the Markov chain 
with the following transition probabilities: the word W 
transforms into W a, where a is one of the generators 
from the set {01, ...,a n } taken with the uniform prob- 
ability 11. For any particular T-step random walk on 
the semigroup B+ we construct the word (the matrix) 
Wt{B^ +1 } and compute the associated largest eigen- 
value At(u) which in the terminology of the work Q 
is called the braiding factor. Then we extract the max- 
imal Lyapunov (" braiding" ) exponent by taking the 
limit 



At = lim 



In \A T (u) 
In u 



Note the difference between Q and the definition of 
the braiding factor in Q. Actually, to find the maximal 
Lyapunov exponent we take the limit with respect to the 
dummy variable u which enters in the definition of the 
matrix representation of the braid group generator (|2|). 
while in [6| the limit is taken with respect to the " time" 
i.e. the length of the word, T. Below we present the argu- 
ments which justify our definition making it topologically 
and physically transparent. 

Since the word Wt is random, the associated Lyapunov 
exponent At is also a random value. Now we proceed as 
it has been stated above. We average A{Wr} over the 
ensemble of all possible sequences Wt and find the mean 
value, A(T). Then we compute the variance, u 2 (T) = 

(A{Wt} — A(T)) 2 and consider its dependence against 
T. The function u(T) is the central result of our letter. 

The matrix representation (|2|) introduced above per- 
mits us to link the topological problems with the " Tetris- 
like" models models of ballistic deposition. Here the no- 
tion of the random walk on the semigroup B^ becomes 
very useful. On one hand, we see from Figjl] that the 
generator <7j has the clear topological meaning, thus the 
information about the degree of entanglement of threads 
(i.e. about the braiding factor) is encoded in the word 
Wt- On the other hand, the word written in terms of 
generators of B+ can be interpreted as a heap (or pile) 
obtained by a 1+1 sequential ballistic deposition process 
in a bounding box. The matrix representation © al- 
lows to determine the corresponding height profile. Let 
us discuss this connection in more details. 
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A standard one-dimensional ballistic deposition model 
with next-ncarcst-neighboring (NNN) interactions is for- 
mulated as follows (for more details, see Refs.(l^4l6j). 
Consider a box divided in n columns (of unit width each) 
enumerated by an index i (i = 1, 2, n). The free con- 
ditions are assumed for left and right boundaries. At the 
initial time moment, (t = 0), the system is empty. Then, 
at each tick of the clock, t = 1,2,...,T, we deposit an 
elementary cell ("particle") of unit height and width in 
a randomly chosen column, i. Suppose that the distribu- 
tion on the set of columns is uniform. Define the height, 
hi(t), in the column i at time moment t. Assume now, 
as it is depicted in Fig[2J that the cells in the nearest- 
neighboring columns interact in such a way that they 
can only touch each other by corners, but never by their 
vertical sides. This implies that after having deposited 
a particle to the column i, the height of this column is 
modified according to the following rule: 

tn(t + 1) = max[hi-i(t), hi(t), h i+ i(t)] + 1 (5) 

If at the time moment t nothing is added to the column 
i, its height remains unchanged: hi(t + 1) = hi(t). A set 
of deposited particles forms a heap (a pile) as shown in 
FigH 
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Figure 2: The heap created by sequential deposition of blocks. 
The interactions of two sequential blocks (1st and 2nd) are 
shown. 

The height hi(t) (1 < i < n) in the column i at time 
t can be easily computed using the matrix representa- 
tion similar to The deposition event in the column 
i means the application of the generator gi with the fol- 
lowing matrix representation: 
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Graphically the deposition process is associated with 
a sequence of "dropping events" of elementary blocks- 
generators gi (1 < i < n). Since we rise a heap by se- 
quential multiplication of matrices ((HJ), we arrive at the 
matrix Ut(ti,u) ='Ilt=i.9it : (1 — ^ — T). Each matrix 
element of Ut(ti,u) is a polynomial of the variable u. 
Take now a vector column a(t = 0) = (oi, ...,a„) where 
the components of a(t = 0) arc distinct nonzero values. 
The local heights h(T) = (/ii(T), K{T)) we extract 
as follows: 



h(T) 



lim 

u— >oo 



ln[U T (n,u) a(i = 0)] 
In it 



(7) 



Let us demonstrate that the definition of heights in ©- 
([7]) is consistent with the updating rules ([3]) for the NNN- 
ballistic deposition process. Using (|6|) we can write the 
recursion relation for the dynamics of the vector a(i): 

a(t + 1) = giMt>, («t6[l,n]) (8) 

If it = i, then 

di(t + 1) = ua,i-\(t) + udiit) + ua,i + i{t) (9) 

Supposing that cii(t) = u hi ^ and substituting this ansatz 
into ©, we get: 



hi(t+l) 



= (7ii_i(t) + l)lnu , Jhi(t) + l)\nu , Jh i+1 (t) + l) In u 



In u 



(10) 



In the limit u — > oo Ea.pU]) coincides with Eq.©. 

The similarity between the representations (J2J and (|6]) 
allows us to expect that the sequential multiplication of 
braid semigroup generators can also be interpreted a sort 
of ballistic deposition process. To see that let us first 
make the following temporary replacement for the 3x3 
block A in 



A = 
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(11) 



By this exchange we ensure for A the absence of any 
cancellations of terms of polynomials in u in course of 
matrix multiplications since all matrix elements stay pos- 
itive. We show later that even for the true block A such 
cancellations are exponentially rare do not change the 
obtained results. Now the multiplication of matrices (0) 
with the replacement A — ^ A has the meaning of the 
heap construction shown in Fig|3] with slightly modified 
local interactions. One can imagine that every elemen- 
tary block has "sticky" right-hand side as it is depicted 
in FigEl 

The corresponding height Hi (t) after depositing a par- 
ticle to the column i at time t is modified as follows: 



(6) H t {t + 1) = maxfffj-iW, H^t) + 1, H i+1 (t) + 1} (12) 
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Figure 3: The heap created by blocks associated with gener- 
ators of the semigroup B^ +1 . The "one-sided-sticky" inter- 
actions of two sequential blocks (1st and 2nd) are shown. 



sign in the block A in . The possible cancellations are 
due to the commcnsurability effects and are negligible for 
generic initial set a(t = 0) = (ai, •■•,a n ) (for instance, if 
o-i, 1 < i < n, are distinct real values from the interval 
]0,1[). The corresponding scaling dependence is shown 
in FigH] (curves 5 and 6). 




log (tln u ) 



If at the time moment t nothing is added to the column 
i, its height remains unchanged: Hi(t + 1) = Hi(t). 

Our first result concerns the computation of the height 
profiles for deposition processes depicted in FigJ5] and 
Fig El The main quantity of interest is the variance 
Wh(T,n) defined for the standard NNN-deposition pro- 
cess as follows: 



w h {T,n) 



1 



,1/2 



n 

E 

.i=i 



(h-hi)' 
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(13) 

The mean value H and the variance wh {T, n) for the bal- 
listic deposition of one-sided-sticky particles (see FigJ3]) 
are defined in the same way. We plot in Fig@] the de- 
pendence wn(t) in the coordinates (w/n 1 / 2 , t/n 3 / 2 ) im- 
plying the known KPZ scaling wh = n 1 / 2 f(t/n 3 / 2 ) with 
f(u) ~ m 1 / 3 for u <C 1 and f(u) ~ const for u^$> 1. 

Comparing figures FigJ5]and Fig(3]we can expect that 
the modifications of local NNN-interactions would not 
affect the scaling behavior of the width of the height dis- 
tribution in a growing heap defined by (1121) . This is ac- 
tually so and as it can be seen by comparing curves (1,2) 
with (3,4) in FigJU Thus, we recover the KPZ scaling for 
the variance WH(t,n) (as the function of "time" t. 

Now we return to the random multiplication of ma- 
trices ([2]) - generators of the braid semigroup, Wt —■ 
ritli a n '■ and numerically compute the associated 



"height", Y=( yi (T) 
Y(T) 



lim 

l— f oo 



y n (T)), defined similarly to (0 
ln[W T {n,u)a(t = 0)] 



(14) 



The interpretation of Y as a hight is still valid however 
some care should be taken to check the absence of statis- 
tically significant amount of cancellations due to the " — " 



Figure 4: KPZ scaling in: the standard NNN-ballistic depo- 
sition (curves 1 and 2); the ballistic deposition of "one-sided- 
sticky" particles (curves 3 and 4); multiplication of braid 
semigroup generators (curves 5 and 6). 

The obtained KPZ scaling for the "heights" depicted 
in FigH] associated with the random multiplication of 
braiding matrices indicates the possibility of similar be- 
havior for the braiding exponent At defined in Q. To 
check that we have multiplied sequentially T randomly 
taken (with uniform distribution // = ^rj) matrices — 
generators of the braid semigroup — from the set 
{<7i, cr n -i} and obtained the associated "word" in the 
matrix representation Wt- Then we have computed nu- 
merically the highest eigenvalue A{WV} and extracted 
the braiding exponent by taking a limit u — > oo as indi- 
cated in (j4]). The results are summarized in FigJS] where 
the KPZ scaling is clearly seen. 
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Figure 5: KPZ scaling in double logarithmic scale for the 
variances in the multiplication of matrices @ with the "true" 
block A. 
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Few words should be said about the procedure of nu- 
merical computation of the resulting matrix Wt and the 
extraction of the highest eigenvalue A{Wr}- The matrix 
multiplication is performed exactly in a symbolic form 
and each matrix element is the polynomial of some de- 
gree of u. Then we have left in each matrix element, Wij , 
only the monomial of highest degree of u. Since later we 
have taken the limit u — > oo onl y su ch monomials dom- 
inate. By Gershgorin theorem |T7| we have found the 
disk of radius r < r(j) = J^Li where the highest 

eigenvalue is located and we have localized the position 
of this highest eigenvalue. We have also checked each 
matrix Wt on the existence of zero's rows. No such rows 
have been detected at least for T < 2 x 10 6 and n < 400. 

The model described here can be extended to the case 
of a " symmetric random stirring protocol" , when the 



symmetric random walk on the full braid group B n+ \ is 
considered. For this situation we expect the same KPZ 
scaling for the variance of the braiding exponent, how- 
ever the description of this process is much more involved 
because the random walk on the full braid group is not a 
Markovian process. Namely, the cancellation in the word 
two sequentially added opposite generators at times t and 
t + 1, like Ui (t) err 1 (t + 1 ) = e demands the knowledge of 
the full prehistory since the initial time moment t = 0. 
This makes the exact analytic approaches to the random 
walks on the full braid group rather doubtful. Neverthe- 
less, we hope that some approximative methods like the 
one developed in [l8[ should still be valid and can be used 
for the investigation of random mixing. 

The authors are grateful to S. Majumdar for valuable 
discussions. 
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